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ABSTRACT 

We obtain new age and mass estimates for the star clusters in M82’s fossil starburst 
region B, based on improved fitting methods. Our new age estimates confirm the peak 
in the age histogram attributed to the last tidal encounter with M81; we find a peak 
formation epoch at slightly older ages than previously published, log(tpeak/yr) = 9.04, 
with a Gaussian a of Alog(twidth) = 0.273. The actual duration of the burst of clus¬ 
ter formation may have been shorter because uncertainties in the age determinations 
may have broadened the peak. Our improved mass estimates conhrm that the (initial) 
masses of the M82 B clusters with V < 22.5 mag are mostly in the range lO'^ — 10 ®Mq, 
with a median mass of Md = 1.08 x 10®M©. The formation history of the observed 
clusters shows a steady decrease towards older ages. This indicates that cluster dis¬ 
ruption has removed a large fraction of the older clusters. 

Adopting the expression for the cluster disruption time-scale of tdis(Af) = 
t 4 *®(M/lO^M 0 )'>' with 7 ~ 0.62 (Paper I), we find that the ratios between the real 
cluster formation rates in the pre-burst phase (log(t/yr) > 9.4), the burst-phase 
(8.4 < log(t/yr) < 9.4) and the post-burst phase (log(t/yr) < 8.4) are about 1:2:^. 
The formation rate during the burst may have been higher if the actual duration of 
the burst was shorter than adopted. 

The mass distribution of the clusters formed during the burst shows a turnover at 
log(Mci/M 0 ) ~ 5.3 which is not caused by selection effects. This distribution can be 
explained by cluster formation with an initial power-law mass function of slope a = 2 
up to a maximum cluster mass of Mmax = 3 x 1 O®M 0 , and cluster disruption with a 
normalisation time-scale f^’^/tburst = (3.0 ±0.3) x 10“^. For a burst age of 1 x 10® yr, 
we find that the disruption time-scale of a cluster of 10^ Mq is ^ 3 x 10^ years, 
with an uncertainty of approximately a factor of two. This is the shortest disruption 
time-scale known in any galaxy. 

Key words: diffusion - galaxies: individual: M82 - galaxies: starburst - galaxies: 
star clusters 


1 INTRODUCTION 

1.1 Multiple starbursts in M82 

M82 is often regarded as the “prototype” starburst galaxy. 
Observations over the entire wavelength range, from radio 
waves to X-rays (see, e.g., Telesco 1988 and Rieke et al. 
1993 for reviews), seem to support a scenario in which tidal 
interactions, predominantly with its large neighbour M81, 
triggered intense star formation in the centre of this small 
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irregular galaxy during the last several 100 Myr. The result¬ 
ing starburst, with the high star formation rate of ~ IOMq 
yr“^, has continued up to about 50 Myr ago (e.g., O’Connell 
& Mangano 1978, Rieke et al. 1993). Energy and gas ejec¬ 
tion from supernovae and combined stellar winds drive a 
large-scale galactic wind along the minor axis of M82 (e.g., 
Lynds & Sandage 1963, McCarthy et al. 1987, Shopbell & 
Bland-Hawthorn 1998). 

All of the bright radio and infrared sources associated 
with the active starburst are found in a small region within a 
radius of ~ 250 pc from the galaxy’s centre, but most of this 
volume is heavily obscured by dust at optical wavelengths. 
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However, there is now compelling evidence that the ac¬ 
tive starburst was not the only major starburst event to 
have occurred in M82. A region at ^ 400— 1000 pc from the 
centre, “M82 B” (nomenclatnre from O’Connell & Mangano 
1978), has the high surface brightness and spectral features 
expected for an ancient starburst with an age in excess of 
several 100 Myr and an amplitude similar to the active burst 
(O’Connell & Mangano 1978, Marcum & O’Connell 1996, de 
Grijs et al. 2001, hereafter dCOG). Its spectral features re¬ 
semble the characteristics of the anomalous “E-l-A” spectra 
found in distant galaxy clusters, which are generally inter¬ 
preted as the signature of a truncated bnrst of star forma¬ 
tion that occurred 100 — 1000 Myr earlier (e.g., Couch & 
Sharpies 1987, Dressier & Gunn 1990, Couch et al. 1998). 
It is thought that these starbursts are closely related to the 
process by which disc galaxies are converted into elliptical or 
lenticular galaxies and they are thought to result from tidal 
interactions, mergers, or perhaps ram-pressure stripping by 
the intergalactic medium (Butcher & Oemler 1978, Oemler 
1992, Barger et al. 1996). 

The significance of the detailed study of M82’s star- 
burst environment lies, therefore, in the broader context of 
galaxy evolution. Starbursts of this scale are likely to be 
common features of early galaxy evolution, and M82 is the 
nearest analogue to the sample of star-forming galaxies re¬ 
cently identified at high redshifts {z > 3; Steidel et al. 1996, 
Giavalisco et al. 1997, Lowenthal et al. 1997). In addition, 
M82 affords a close-up view not only of an active starburst 
- in M82 A, G and E - but also, in region B and elsewhere, 
of the multiple post-burst phases. 

In a recent study focusing on the fossil starburst site, 
M82 B, we found a large population of (~ 110) evolved com¬ 
pact star clusters (dGOG), whose properties appear to be 
consistent with them being evolved (and therefore faded) 
counterparts of the young (super) star clusters detected in 
the galaxy’s active core (O’Gonnell et al. 1995). Based on 
broad-band optical and near-infrared colours and compari¬ 
son with stellar evolutionary synthesis models, we estimated 
ages for the M82 B cluster population from ~ 30 Myr to over 
10 Gyr, with a peak near 650 Myr (see also de Grijs 2001). 
About 22 per cent of the clusters in M82 B are older than 
2 Gyr, with a roughly flat distribution to over 10 Gyr. Very 
few clusters are younger than 300 Myr. The full-width of 
the peak at 650 Myr is ~ 500 Myr, but this is undoubt¬ 
edly broadened by the various uncertainties entering the 
age-dating process. The selection bias of the star clusters 
in dGOG is such that the truncation of cluster formation 
for ages < 300 Myr is better established than the roughly 
constant formation rate at ages > 2 Gyr. 

Thus, we suggested steady, continuing cluster formation 
in M82 B at a very modest rate at early times (> 2 Gyr ago) 
followed by a concentrated formation episode lasting from 
400-1000 Myr ago and a subsequent suppression of cluster 
formation (dGOG). 

1.2 Cluster disruption time-scales 

In the determination of the cluster formation history from 
the age distribution of magnitude-limited cluster samples, 
cluster disruption must be taken into account. This is be- 
canse the observed age distribution is that of the surviving 
clusters only. Therefore, in this study of the evolution of 


the star cluster system in M82’s fossil starburst region, we 
determine both its cluster formation history and the char¬ 
acteristic disruption time-scales. 

The dynamical evolution of star clusters is determined 
by a combination of internal and external time-scales. The 
free-fall and two-body relaxation time-scales, which depend 
explicitly on the initial cluster mass density (e.g., Spitzer 
1957, Ghernoff & Weinberg 1990, de la Fuente Marcos 1997, 
Portegies Zwart et al. 2001), affect the cluster-internal pro¬ 
cesses of star formation and mass redistribution through en¬ 
ergy equipartition, leading to mass segregation and, even- 
tnally, core collapse (see, e.g., de Grijs et al. [2002a,b] for 
a detailed description of mass segregation effects in young 
star clusters). While the internal relaxation process will, over 
time, eject both high-mass stars from the core (e.g., due to 
interactions with hard binaries; see Brandi et al. 2001, de 
Grijs et al. 2002a) and lose lower-mass stars from its halo 
through diffusion (e.g., due to Roche-lobe overflow), the ex¬ 
ternal process of tidal disruption and stripping by the sur- 
ronnding galactic field is in general more important for the 
discussion of the disruption of star clusters. 

Tidal disruptive processes are enhanced by “normal” 
stellar evolutionary effects such as mass loss by winds and/or 
supernova explosions, which will further reduce the stellar 
density in a cluster, and thus make it more sensitive to ex¬ 
ternal tidal forces. 

From the bimodal age distribution of (young) open and 
(old) globular clusters in the Milky Way, Oort (1957) con¬ 
cluded that disruption of Galactic star clusters must occur 
on time-scales of ~ 5 x 10® yr. Around the same time, Spitzer 
(1957) derived an expression for the disruption time scale 
as a function of a cluster’s mean density, pc (M© pc“®): 
tdis = 1.9 X 10®Pc yr, for 2.2 < pc < 22 Mq pc“®. More 
advanced recent studies, based on A-body modeling, have 
shown that the cluster disruption time-scale is sensitive to 
the cluster mass, the fraction of binary (or multiple) stars, 
and the initial mass function (IMF) adopted (e.g., Ghernoff 
& Weinberg 1990, de la Fuente Marcos 1997). 

Boutloukos & Lamers (2001, 2002) derived an empiri¬ 
cal relation between the disruption time and a cluster’s ini¬ 
tial mass for the Milky Way, the Small Magellanic Cloud 
(SMC), M33 and the inner spiral arms of M51. They showed, 
based on an analysis of the mass and age distributions of 
magnitude-limited samples of clusters, that the empirical 
disruption time of clusters depends on their initial mass (Mi) 
as 

tdis = 4“(Mi/10‘‘M©)^ (1) 

where 14 '° is the disruption time of a cluster of initial mass 
Mi = IGMq. The value of 7 is approximately the same 
in these four galaxies, 7 = 0.62 ± 0.06. However, the char¬ 
acteristic disruption time-scale 14 ’” is widely different in the 
different galaxies. The disruption time-scale is longest in the 
SMC (~ 8 Gyr) and shortest in the inner spiral arms of M51 
(^ 40 Myr). We will derive the disruption time-scale of clus¬ 
ters in M82 B and allow for disruption in the determination 
of the cluster formation history. 

In Section 2, we will summarise the observations on 
which our discussion is based, and in Section 3 we ob¬ 
tain new age and mass estimates for the clusters in M82 
B, based on the 3/2DEF method (3/2-dimensional energy¬ 
fitting method; Bik et al. 2002) of fitting the observed spec- 
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tral energy distribution to that of cluster evolution models. 
In Section 4 we derive the cluster formation history and the 
characteristic cluster disruption time-scale. In Section 5 we 
discuss the mass distribution of the clusters formed during 
the burst. The methods and assumptions are discussed in 
Section 6; finally, we will summarise and conclude the paper 
in Section 7. 


2 OBSERVATIONS 

The observations on which our analysis of the M82 clus¬ 
ter formation history is based were described in detail in 
dGOG. Briefly summarised, we observed region B with the 
Wide Field Planetary Camera 2 (WFPC2) on board the 
Hubble Space Telescope (HST) through the F439W, F555W 
and F814W filters, roughly corresponding to the standard 
Johnson-Cousins B, V and / passbands, respectively. We 
obtained WFPC2 observations using two pointings, so that 
the entire region B was covered by the Planetary Camera 
(PC) chip, the highest-resolution optical imaging instru¬ 
ment available on board HST at that time, with a pixel 
size of 0.0455”. The effective integration times used for the 
F439W, F555W, and F814W observations were 4400s, 2500s 
and 2200s, respectively, for the western half of M82 B and 
4100s, 3100s and 2200s, respectively, for the eastern half 
closest to the galaxy’s starburst core. 

In addition, we imaged the entire region with HSTs 
Near-Infrared Camera and Multi-Object Spectrometer 
(NICMOS) (Camera 2; pixel size 0.075”) in both the FllOW 
and F160W filters (comparable to the Bessell Jand //filters, 
respectively), in a tiled pattern of 2 x 4 partially overlapping 
exposures. The integrations, with effective integration times 
of 768s for each of the eastern and western halts of M82 
B and each filter, were taken in MULTIACCUM mode to 
preserve dynamic range and to correct for cosmic rays. 

We subsequently obtained integrated photometry for 
the extended objects (i.e., star clusters) in M82 B down to 
a 50 per cent completeness limit of V ~ 23.3 mag. Because 
of the highly variable background and the numerous spuri¬ 
ous features due to dust lanes and background variations, 
we decided to select only genuine star cluster brighter than 
V — 22.5 for further analysis, corresponding to close to 100 
per cent completeness (dGOG). 


3 OBSERVED CLUSTER MASS AND AGE 

DISTRIBUTIONS 

For the final 113 objects obtained following the procedures 
outlined in the previous section, dGOG used (B — V) versus 
(V — I) colour-colour diagrams to disentangle age and ex¬ 
tinction effects, since the age and extinction vectors are not 
entirely degenerate for this choice of optical colours. 

Accurate age and mass determinations are essential for 
the analysis of cluster formation and disruption time-scales 
performed in this paper. Combining the observed luminosi¬ 
ties of the M82 B star cluster population with the appropri¬ 
ate, age-dependent mass-to-light ratios from spectral syn¬ 
thesis models provided us with photometric mass estimates 
with an accuracy of well within an order of magnitude. In¬ 
dependent dynamical mass estimates from high-resolution 


spectroscopy are available only for a few of the most lumi¬ 
nous super star clusters seen in the nearest starburst galax¬ 
ies (M82, the Antennae, NGC 1569 and NGC 1705), and are 
approximately 10®M© (Ho & Filippenko 1996a,b, Smith & 
Gallagher 2000, Mengel et al. 2002). Because of the prox¬ 
imity of M82, we have been able to probe the young cluster 
population in M82 B to fainter absolute magnitudes, and 
thus lower masses, than has been possible before in other 
starburst galaxies. Other young star cluster samples are bi¬ 
ased towards high masses by selection effects due to their 
host galaxies’ greater distances. 

Keeping the importance of the age dependence on our 
photometric mass estimates in mind, we decided to redeter¬ 
mine the M82 B cluster ages using the full parameter space 
available. We obtained more accurate age determinations by 
also including the NICMOS observations, and solving for the 
best-fitting ages matching the full spectral energy distribu¬ 
tions (SEDs) of our sample clusters, from the F439W to the 
F160W passband. We compared the observed cluster SEDs 
with the model predictions for an instantaneous burst of star 
formation, assuming a Salpeter IMF from 0.1 — lOOM© with 
power-law slope oimf = 2.35 to obtain our new estimates for 
the cluster age t, initial mass M^i and extinction E{B — V); 
for the latter we adopted the (Galactic) extinction law of 
Scuderi et al. (1996). For ages t < 10® yr we used the Star- 
burst99 models (Leitherer et al. 1999), while for older ages 
we used the most recent single stellar population models by 
Bruzual & Chariot (2000, hereafter BCOO). Examples of our 
model fits to the observed cluster SEDs are shown in Fig. 1. 

We realise that recent determinations of the stellar 
IMF deviate significantly from a Salpeter-type IMF at low 
masses, in the sense that the low-mass stellar IMF is sig¬ 
nificantly flatter than the Salpeter slope. The implication 
of using a Salpeter-type IMF for our cluster mass deter¬ 
minations is therefore that we have overestimated the indi¬ 
vidual cluster masses (although the relative mass distribu¬ 
tion of our entire cluster sample remains unaffected). There¬ 
fore, we used the more modern IMF of Kroupa, Tout & 
Gilmore (1993, hereafter KTG) to determine the correc¬ 
tion factor, C, between our masses and the more realistic 
masses obtained from the KTG IMF (both normalised at 
l.OM©). This IMF is characterised by slopes of a = —2.7 
for m > l.OM©, a = —2.2 for 0.5 < m/MQ < 1.0, and 
— 1.85 < a < —0.70 for 0.08 < m/MQ < 0.5. Depending on 
the adopted slope for the lowest mass range, we have there¬ 
fore overestimated our individual cluster masses by a factor 
of 1.70 < C < 3.46 for an IMF containing stellar masses in 
the range 0.1 < m/MQ < 100 . 

The fitting of the observed cluster SEDs to the 
Starburst99 and BCOO models was done using a three- 
dimensional maximum likelihood method, 3/2DEF, with the 
initial mass Mi, age and extinction E{B — V) as free param¬ 
eters (see Bik et al. 2002). For the 46 clusters with upper 
limits in one or more filters we used a two-dimensional maxi¬ 
mum likelihood fit, using the extinction probability distribu¬ 
tion for E{B — V). This distribution was derived for the clus¬ 
ters with well-defined SEDs over the full wavelength range 
(see Bik et al. 2002). We obtained reliable age estimates 
based on the full SED modeling for 81 of the 113 clusters. 
The new age and mass estimates for these 81 clusters are 
listed in Table 1, where the individual clusters are identified 
by their dGOG ID. The dGOG mass estimates in Table 1 


© 2001 RAS, MNRAS 000, 000-000 



4 R. de Grijs, N. Bastian, and H.J.G.L.M. Lamers 


Table 1. 


Notes: “ 
cluster 


New age, mass and extinction determinations for the M82 B clusters, and comparison with previously published values. 


This paper dGOG 


Object 

min 

log(Age/yr) 

best 

max 

log(m/MQ)“ 

min best max 

E(S 

min 

- V) (mag) 
best max 

Iog(Age/yr) 

log{m/MQ) 

E(B - V) 
(mag)*’ 

Bl-01 

9.16 

9.88 

10.05 

5.20 

5.62 

5.74 

0.00 

0.02 

0.20 

10.24 

5.92 

0.31 

Bl-02 

8.96 

9.28 

9.36 

4.66 

4.90 

4.98 

0.00 

0.00 

0.16 

9.27 

4.90 

0.04 

Bl-04 

5.10 

8.46 

8.81 

4.14 

5.21 

5.34 

0.06 

0.18 

0.70 

8.78 

5.25 

0.58 

Bl-05 

9.74 

10.05 

10.15 

5.76 

5.95 

6.03 

0.00 

0.00 

0.06 

10.14 

6.04 

0.18 

Bl-06 

5.10 

6.58 

8.76 

4.33 

4.53 

5.54 

0.00 

0.52 

0.60 

8.67 

5.46 

0.44 

Bl-07 

8.96 

9.01 

9.16 

5.53 

5.60 

5.69 

0.00 

0.06 

0.08 

9.03 

5.53 

0.12 

Bl-08 

8.51 

8.66 

8.66 

5.14 

5.24 

5.24 

0.00 

0.02 

0.04 

8.68 

5.24 

0.09 

Bl-09 

5.10 

8.91 

9.01 

3.14 

4.24 

4.34 

0.00 

0.00 

0.72 

8.89 

4.21 

0.00 

Bl-10 

5.10 

8.71 

8.86 

3.14 

4.33 

4.41 

0.00 

0.04 

0.64 

8.76 

4.33 

0.18 

Bl-11 

8.71 

8.86 

8.86 

5.36 

5.38 

5.41 

0.00 

0.00 

0.08 

8.84 

5.36 

0.09 

Bl-12 

8.71 

9.01 

9.21 

4.81 

4.88 

5.02 

0.00 

0.06 

0.22 

9.03 

4.81 

0.06 

Bl-14 

8.51 

8.86 

9.01 

5.11 

5.24 

5.28 

0.00 

0.10 

0.22 

8.89 

5.13 

0.34 

Bl-16 

8.56 

9.95 

10.30 

4.94 

5.63 

5.87 

0.00 

0.04 

0.50 

10.00 

5.63 

0.40 

Bl-17 

8.66 

9.01 

9.16 

4.47 

4.56 

4.68 

0.00 

0.04 

0.24 

9.03 

4.50 

0.17 

Bl-18 

8.86 

8.91 

8.91 

5.54 

5.56 

5.56 

0.00 

0.00 

0.02 

8.85 

5.50 

0.15 

Bl-20 

9.01 

9.16 

9.16 

5.48 

5.65 

5.65 

0.00 

0.00 

0.08 

9.10 

5.60 

0.07 

Bl-22 

8.56 

8.76 

8.86 

5.23 

5.29 

5.33 

0.00 

0.04 

0.12 

8.78 

5.26 

0.25 

Bl-24 

8.66 

8.86 

8.91 

5.15 

5.20 

5.23 

0.00 

0.02 

0.12 

8.83 

5.15 

0.13 

Bl-25 

8.91 

9.01 

9.26 

5.00 

5.10 

5.23 

0.00 

0.10 

0.16 

9.05 

4.97 

0.16 

Bl-26 

5.10 

8.61 

8.66 

3.09 

4.26 

4.32 

0.00 

0.00 

0.56 

8.68 

4.33 

0.00 

Bl-27 

8.76 

8.86 

8.91 

4.44 

4.49 

4.51 

0.00 

0.02 

0.08 

8.86 

4.46 

0.00 

Bl-28 

8.86 

8.86 

8.86 

5.97 

5.97 

5.97 

0.04 

0.04 

0.04 

8.86 

5.90 

0.13 

Bl-29 

5.10 

8.76 

8.91 

3.52 

4.71 

4.78 

0.00 

0.06 

0.68 

8.81 

4.68 

0.25 

Bl-32 

8.76 

9.01 

9.01 

4.64 

4.69 

4.72 

0.00 

0.02 

0.14 

8.98 

4.65 

0.10 

Bl-34 

5.10 

6.84 

8.36 

3.31 

3.39 

4.44 

0.00 

0.12 

0.46 

7.55 

3.94 

0.35 

Bl-35 

8.61 

8.86 

9.01 

5.02 

5.08 

5.12 

0.00 

0.08 

0.20 

8.88 

4.98 

0.24 

Bl-36 

8.96 

10.14 

10.30 

5.13 

5.90 

6.01 

0.00 

0.00 

0.34 

10.14 

5.91 

0.24 

Bl-37 

9.16 

9.40 

9.44 

5.75 

5.88 

5.91 

0.00 

0.00 

0.10 

9.32 

5.80 

0.18 

Bl-38 

5.10 

8.81 

8.96 

3.57 

4.63 

4.71 

0.00 

0.04 

0.68 

8.83 

4.60 

0.16 

Bl-39 

8.81 

9.34 

9.44 

4.29 

4.58 

4.70 

0.00 

0.00 

0.28 

9.28 

4.53 

0.15 

Bl-40 

5.10 

8.66 

8.76 

2.79 

3.97 

4.03 

0.00 

0.02 

0.60 

8.68 

3.97 

0.13 

Bl-43 

9.01 

9.48 

9.68 

4.35 

4.69 

4.87 

0.00 

0.00 

0.20 

9.70 

4.90 

0.11 

B2-01 

9.01 

9.94 

10.25 

4.43 

5.09 

5.34 

0.00 

0.00 

0.28 

10.00 

5.16 

0.00 

B2-04 

9.01 

9.23 

9.28 

5.23 

5.44 

5.47 

0.00 

0.02 

0.14 

9.32 

5.51 

0.15 

B2-05 

5.10 

8.46 

8.71 

4.14 

5.21 

5.35 

0.00 

0.08 

0.62 

8.66 

5.28 

0.32 

B2-07 

8.66 

8.96 

9.41 

5.41 

5.55 

5.78 

0.00 

0.18 

0.34 

9.19 

5.56 

0.43 

B2-08 

8.66 

9.01 

10.26 

5.49 

5.60 

6.35 

0.00 

0.32 

0.50 

9.72 

5.89 

0.77 

B2-12 

9.74 

9.76 

9.80 

6.86 

6.88 

6.91 

0.00 

0.00 

0.00 

9.70 

6.83 

0.22 

B2-13 

9.30 

9.83 

9.90 

5.38 

5.71 

5.77 

0.00 

0.00 

0.12 

9.76 

5.66 

0.19 

B2-14 

8.71 

9.44 

10.14 

4.75 

5.13 

5.63 

0.00 

0.08 

0.40 

10.13 

5.64 

0.43 

B2-17 

8.96 

9.48 

9.80 

5.21 

5.56 

5.81 

0.00 

0.02 

0.26 

9.46 

5.52 

0.26 

B2-18 

5.10 

7.49 

8.66 

4.09 

4.79 

5.44 

0.08 

0.26 

0.74 

9.23 

5.69 

1.03 

B2-21 

6.86 

7.00 

7.60 

4.46 

4.57 

5.17 

0.00 

0.04 

0.24 

8.90 

5.82 

1.13 

B2-22 

8.96 

10.20 

10.30 

4.53 

5.36 

5.47 

0.00 

0.00 

0.36 

10.26 

5.39 

0.15 

B2-25 

5.10 

8.51 

8.81 

3.90 

5.01 

5.14 

0.00 

0.10 

0.64 

8.71 

5.05 

0.29 

B2-26 

8.71 

8.76 

8.76 

5.60 

5.62 

5.62 

0.02 

0.02 

0.04 

8.77 

5.63 

0.14 

B2-28 

5.10 

7.49 

8.36 

3.51 

4.09 

4.65 

0.00 

0.00 

0.48 

7.26 

3.91 

0.60 
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based on a Salpeter-type IMF; comparison with results from more modern IMFs (e.g., KTG) implies that we have overestimated our individual 
masses by a factor of about 1.7—3.5, depending on the adopted IMF slope for the lowest stellar masses (see Sect. 3); ^ E(S — V) = Ay /3.I. 
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Figure 1. Examples of a few fits of the Starburst99 and BCOO models (dashed lines) to the observed cluster SEDs sampled in the 
B, V, I, J, and H passbands (data points with error bars), representative of clusters throughout the entire age range (age increases from 
panel a to f). The best-fitting model SEDs are based on the 3/2DEF maximum likelihood method. 

were determined by combining the absolute magnitudes of justified. We should caution, however, that we have only 

the clusters with their age-dependent mass-to-light ratios. been able to sample the star clusters located close to the 

surface of region B, as evidenced by the derived extinction 
Our analysis in both dGOG and the present paper is estimates: dGOG find Av < 1 mag in general, while our 

based on the assumption of solar metallicity, which should new determinations (this paper) restrict the extinction even 

be a reasonable match to the young objects in M82 (Gal- more, to Av < 0.2 mag for the clusters with well-determined 

lagher & Smith 1999; see also Fritze-v. Alvensleben & Ger- ages, and possibly up to Ay ^ 0.40 — 0.55 for some of the 

hard 1994), but the effects of varying the metallicity between others, 

one-fifth and 2.5 times solar are small compared to the pho¬ 
tometric uncertainties (see also dGOG). We therefore assert In Fig. 2 we compare the age and mass estimates from 

that our results are not significantly affected by possibly the more sophisticated 3/2DEF method used in this pa- 

varying metallicities. Ongoing analysis (Parmentier, de Grijs per with the corresponding parameters obtained by dGOG 

& Gilmore, 2002) suggests that this assumption was indeed based on their location in the {B — V) vs. {V — I) di- 
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agrams. The black dots represent clusters for which the 
range between the minimum and maximum plausible ages, 
log(Age[max]) —log(Age[min]) < 1.0; the open circles are ob¬ 
jects with more uncertain age determinations. We conclude 
that the new age and mass determinations are consistent 
with the dGOG values; this consistency is better for the 
photometric cluster mass estimates than for their ages, as 
shown by the larger scatter in Fig. 2a with respect to panel 
b. However, we also note that the age and mass distributions 
obtained for the subsample with well-determined ages and 
for the full sample are internally consistent. 

Fig. 3 shows the distribution of the M82 B clusters in 
the age vs. mass plane. It is immediately clear that the lower 
mass limit increases with increasing cluster age. The vari¬ 
ous (solid, dashed and dotted) lines overplotted on the figure 
show that this is indeed the expected effect of normal evolu¬ 
tionary fading of a synthetic single stellar population of an 
instantaneously formed cluster. For ages < 10® yr, we show 
the unreddened fading line for clusters with a limiting mag¬ 
nitude oiV = 22.5 at the distance of M82 {m—M — 27.8; see 
dGOG), and various choices for the IMF (indicated in the 
figure are the IMF slope oimf, and the lower and upper mass 
cut-offs), predicted by the Starburst99 models (“SB99”). For 
older ages (t > 10® yr), we show its extension predicted by 
the BGOO models. These predicted lower limits agree well 
with our data points. 


4 THE DERIVED CLUSTER FORMATION 

HISTORY 

4.1 Cluster age and mass distributions 

The cluster age and mass distributions for the M82 B 
cluster sample are shown in Figs. 4a and b. Our new 
age estimates confirm the peak in the age histogram at¬ 
tributed to the last tidal encounter with M81 by dGOG. 
For the subsample of clusters with well-determined ages 
(shaded histogram), we find a peak formation epoch at 
log(tpeak/yr) = 9.04 (tpeak — 1.10 Gyr), with a Gaussian 
cr of A log(t„idth) = 0.27, corresponding to a FWHM of 
A log(twidth) = 0.64. The corresponding numbers for the 
full sample are log(tpeak/yr) = 8.97 (tpeak — 0.93 Gyr) for 
the peak of cluster formation, and for a, A log(twidth) = 0.35 
(FWHM, Alog(twidth) = 0.82). In dGOG (see also de Grijs 
2001 ), we concluded that there is a strong peak of cluster 
formation at ~ 650 Myr ago, which have formed over a pe¬ 
riod of ^ 500 Myr, but very few clusters are younger than 
log(t/yr) ~ 8.5 {t ~ 300 Myr). Our new age estimates date 
the event triggering the starburst to be slightly older than 
the value of 650 Myr derived by dGOG. Our estimate of the 
duration of the burst should be considered an upper limit, 
because uncertainties in the cluster age determination may 
have broadened the peak in Fig. 4a. 

At the distance from the centre of region B, one would 
expect M82’s differential rotation (cf. Shen & Lo 1995) to 
have caused the starburst area to disperse on these time- 
scales. The reason why the fossil starburst region has re¬ 
mained relatively well constrained is likely found in the com¬ 
plex structure of the disc. It is well-known that the inner 
^ 1 kpc of M82 is dominated by a stellar bar (e.g.. Wills et 
al. 2000) in solid-body rotation. From observations in other 


galaxies, it appears to be a common feature that central bars 
are often surrounded by a ring-like structure. If this is also 
true for M82, it is reasonable to assume that stars in the ring 
are trapped, and therefore cannot move very much in radius 
because of dynamical resonance effects. The phase mixing 
around the ring might be slow enough for a specific part of 
the ring to keep its identity over a sufficient time so as to 
appear like region B (see de Grijs 2001): if the diffusion ve¬ 
locity around the ring is sufficiently small, any specific region 
would remain self-constrained for several rotation periods. 
In addition, since the density in the region is high (see Sect. 
6 .1), simple calculations imply that the area’s self-gravity 
is non-negligible compared to the rotational shear, therefore 
also prohibiting a rapid dispersion of the fossil starburst re¬ 
gion. 

The improved mass estimates also confirm that the (ini¬ 
tial) masses of the young clusters in M82 B with V < 22.5 
mag are mostly in the range lO'^ — 10® M©, with a me¬ 
dian of 10®M© (dGOG). In fact, based on a Gaussian fit, 
we find that the mean mass of our M82 B cluster sample 
is log(Mci/M©) = 5.03 and 4.88 for the subsample with 
well-determined ages and the full sample, respectively, cor¬ 
responding to Mci = 1.1 X 10® and 0.8x 10®M©, respectively. 

4.2 Cluster disruption and the cluster formation 
history 

In Figs. 4c and d we show the formation rate and the mass 
spectrum of the observed clusters in M82 B. The open cir¬ 
cles were derived for the entire sample, whereas the filled 
data points represent only those clusters for which we could 
obtain good age estimates. These distributions depend on 
the cluster formation history and on the cluster disruption 
time-scale of M82 B. 

As shown by Boutloukos & Lamers (2002, hereafter Pa¬ 
per I), with only a few well-justified assumptions, the mass 
and age distributions of a magnitude-limited sample of star 
clusters in a given galaxy can be predicted both accurately 
and robustly, despite the complex physical processes under¬ 
lying the assumptions (for a full discussion see Paper I). If 
all of the following conditions are met: 

(i) the cluster formation rate, = S x M©“, is 

constant; 

(ii) the slope a of the cluster IMF is constant with 
V(Mei)dM « M©“dM; 

(iii) stellar evolution causes clusters to fade as F\ ~ t~‘’, 
as predicted by cluster evolution models; and 

(iv) the cluster disruption time depends on their initial 
mass as tdis = x (Md/IO'^M©)'’', where is the disrup¬ 
tion time-scale of a cluster with M© = lO’^M©. It is well- 
established, however, that the disruption time-scale does not 
only depend on mass, but also on the initial cluster density 
and internal velocity dispersion (e.g., Spitzer 1957, Cher- 
noff & Weinberg 1990, de la Fuente Marcos 1997, Portegies 
Zwart et al. 2001). Following the approach adopted in Pa¬ 
per I, however, we point out that if clusters are approxi¬ 
mately in pressure equilibrium with their environment, we 
can expect the density of all clusters in a limited volume 
of a galaxy to be roughly similar, so that their disruption 
time-scale will predominantly depend on their (initial) mass 
(with the exception of clusters on highly-eccentric orbits). In 
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Figure 2. Comparison between (a) the age and (b) mass estimates of dGOG based on BVI photometry and new estimates using 
additional JH observations and a full-parameter fitting routine as described in the text. The black dots represent clusters for which the 
total age range obtained is (log(Age[max]) — log (Age [min])) < 1.0; the open circles are objects with more uncertain age determinations. 
They show a larger scatter than the black dots, but both distributions are internally consistent. 


the opposite case that the initial cluster density p depends 
on their mass M in a power-law fashion, e.g., p ^ with 
X being the (arbitrary) power-law exponent, the disruption 
time-scale will also depend on mass if tdis ^ (Paper 

I). 

then it can be shown easily that the age distribution of the 
observed cluster population will obey the following approx¬ 
imate power-law behaviours (see Paper I): 

• dA7i/dt oc for young clusters due to fading; 

• dA7i/dt oc for old clusters due to disruption. 

Similarly the mass spectrum of the observed clusters will be 

• dA7i/dMci oc for low-mass clusters due to 

fading; 

• dVci/dMci oc for high-mass clusters due to dis¬ 

ruption. 

So both distributions will show a double power law with 
slopes determined by a, and 7 . The crossing points are 
determined by the cluster formation rate and the character¬ 
istic disruption time-scale of a cluster with an initial mass 
of tfy 

In Paper I we showed that the observed age and mass 
distributions of the star cluster systems in four well-studied 
galaxies indeed show the predicted double power-law be¬ 
haviour. From the analysis of these observed distributions 
Boutloukos & Tamers (2002) showed that the value of 7 is 
approximately constant, (7 = 0.62 ± 0.06), under the very 
different environmental conditions in the four galaxies, but 
the characteristic disruption time-scales differ significantly 
from galaxy to galaxy. 


For the clusters in M82 B the situation is more com¬ 
plex, because the cluster formation rate was certainly not 
constant. In fact the distribution in Fig. 4c suggests that 
we can distinguish three phases, which we wifi discuss sepa¬ 
rately. 

We will assume that the cluster formation rate has been 
constant within each phase, but may differ strongly among 
the phases. We adopt a cluster IMF of slope a = 2.0 (Harris 
& Pudritz 1994; McLaughlin & Pudritz 1996; Elmegreen & 
Efremov 1997; Zhang & Fall 1999; Bik et al. 2002). We will 
also adopt a slope C, = 0.648 for the evolutionary fading of 
clusters in the V band, derived from the Starburst99 clus¬ 
ter models (see also Paper I). For the mass scaling of the 
disruption time-scale we adopt 7 = 0.62 (Paper I). 

(i) The pre-burst phase, at log(t/yr) > 9.4. The age dis¬ 
tribution of this earliest phase is represented by only three 
age-bins of high age at log(t/yr) > 9.4. The predicted slope 
is (1 — a )/7 = —1.61. We have fitted a straight line through 
the last three age-bins with this slope. We see that the pre¬ 
dicted slope matches the observations within the uncertain¬ 
ties. 

(ii) The burst phase, in the approximate time interval 
8.4 < log(t/yr) < 9.4. Its distribution can be fitted with 
a line of the same predicted slope of —1.61. This is also 
shown in Fig. 4c. The predicted slope fits the observations 
very well. We see that this disruption line of the burst-phase 
is higher than in the pre-burst phase by (0.3T0.1) dex. This 
implies that the cluster formation rate during the burst was 
approximately a factor of 2 higher than before. However, 
this factor depends on the assumed duration of the burst. 
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Figure 3. Distribution of the M82 B clusters in the (age vs. mass) plane. The symbol coding is as in Fig. 2. Overplotted for ages up 
to 1.0 Gyr are the expected detection limits in the (age vs. mass) plane, predicted by Starburst99 for a range of IMFs; for older ages, 
we use the BCOO models for a standard IMF, as indicated in the figure legend. These model predictions are based on a detection limit 
of y = 22.5 and (m — M)m 82 = 27.8, assuming no extinction. For a nominal extinction of Ay = 0.2 mag, expected for the clusters 
with well-determined ages, the detection limit is expected to shift to higher masses by A log(A7(,i/Af0) = 0.08, which is well within the 
uncertainties associated with our mass determinations. The features around 10 Myr are caused by the appearance of red supergiants. 


The duration of the burst may have been shorter than the 
value of Alog(t) = 1.0 suggested by the figure, because un¬ 
certainties in the derived cluster ages may have broadened 
the burst peak. If the duration of the burst was a factor 
X shorter than we adopted, then the overall cluster forma¬ 
tion rate during the burst was a factor X higher than we 
estimated. 

(iii) The post-burst phase, at log(t/yr) < 8.4. The age 
distribution of this most recent, post-burst phase extends 
from 6.0 < log(t/yr) < 8.4. Unless the disruption time 
was extremely short with of order a few Myr (which 
is contradicted by the analysis of the mass distribution in 
the burst, see Sect. 5), the age distribution of the youngest 
bins is governed by the fading of the clusters and the 
detection limit. Despite the rather large uncertainties for 
young ages, log(t/yr) < 7.0, we find a best-fitting slope 


C(1 - a) = -0.57 ± 0.05, so that ( = 0.57 db 0.05 (for¬ 
mal uncertainty, not including the large error bars) for a 
mass-function slope a = 2 (dotted line in Fig. 4c). This is, 
within the rather large uncertainties, similar to the value of 
^ ~ 0.65 expected from the StarburstOO and BCOO models. 
Therefore we have assumed a fading line with the predicted 
slope (((1 — a) = —0.65 for the youngest three age bins. The 
age bins in the range of 7.0 < log(t/yr) < 8.4 fall below the 
extrapolated fading line. The difference is about one dex 
at log(t/yr) ~ 8, although the number of observed clusters 
is small. This suggest that the number of clusters for ages 
f ^ 10^ yr is already affected by disruption. We can fit a 
disruption line with the predicted slope of —1.61 through 
the data points, despite their rather large uncertainties. We 
see that the best-fitting line is displaced by ( — 1.6 ±0.1) dex 
compared to the pre-burst disruption line. This shows that 
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Figure 4. (a) and (b) - Age and mass distributions of the M82 B clusters, redetermined in this paper using improved methods. The 
shaded histograms correspond to the clusters with (log(Age[max]) — log(Age[min]) < 1.0; the open histograms represent the entire cluster 
sample, (c) - The cluster formation rate (in number of clusters per Myr) as a function of age. Open circles: full sample; filled circles: 
clusters with well-determined ages, as above. The vertical dashed lines indicate the age range dominated by the burst of cluster formation. 
The dotted line is the least-squares power-law fit to the fading, non-disrupted clusters, for a constant ongoing cluster formation rate. 
The solid line segments are the disruption lines of clusters formed in the pre-burst phase, during the burst, and in the post-burst phase, 
as described in Section 6.1. (d) - Mass spectrum of the M82 B clusters (number of clusters per mass bin); symbol coding as in panel (c). 


the cluster formation rate in the post-burst phase was ap¬ 
proximately a factor of 40 smaller than during the pre-burst 
phase. 

We conclude from this analysis that the cluster forma¬ 
tion rate during the burst was at least a factor of 2 higher 
than during the pre-burst phase (depending on the duration 
of the burst), and that the cluster formation rate was a factor 
of ~ 40 smaller after the burst than before the burst. This 
is not surprising because the intense episode of cluster for¬ 
mation during the burst will have consumed a large fraction 
of the available number of molecular clouds, leaving little 
material for cluster formation at later times. In fact, it is 
likely that starbursts are strongly self-limited, or quenched. 


by supernova-driven outflows, which remove the remaining 
cool gas from the immediate starburst region (e.g., Cheva¬ 
lier & Clegg 1985, Doane & Mathews 1993). The remarkable 
minor-axis wind in M82 is a dramatic example of this pro¬ 
cess. However, the disturbed conditions near an early burst 
may discourage re-ignition at the same site when cool gas 
inflows resume, shifting the location of active star formation 
(e.g., dCOG). 
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5 THE MASS DISTRIBUTION OF CLUSTERS 

FORMED IN THE BURST. 

In this section we analyse the mass distribution of the clus¬ 
ters that were formed during the burst of cluster forma¬ 
tion, roughly defined to have occurred in the age range 
8.4 < log(t/yr) < 9.4 (see the vertical dashed lines in Fig. 
4c). In Fig. 5 we show their mass distribution. Although the 
M82 B cluster system does not contain large numbers of 
very massive clusters, its mass distributions shows a clear 
increase in cluster numbers from the high-mass end towards 
lower masses, for logMd/MQ > 5. For lower masses, how¬ 
ever, the number of clusters decreases rapidly. If the mass 
spectrum of young, newly formed clusters resembles a power- 
law distribution, this turn-over at log Mci/Mq ~ 5 is not ex¬ 
pected, unless disruption effects have preferentially removed 
the lower-mass clusters from our magnitude-limited sample. 

In dGOG we concluded that the distribution of clus¬ 
ter luminosities (and therefore the equivalent mass distri¬ 
bution) in M82 B - corrected to a common age of 50 Myr 
- shows a broader and flatter mass distribution than typ¬ 
ical for young cluster systems, although this distribution 
is subject to strong selection effects. For the small age 
range considered here, the selection limit in observable mass 
at log{M/ Mq) ~ 4.2 imposed by the brightness limit at 
V = 22.5 mag (Fig. 3) occurs at almost an order of mag¬ 
nitude lower masses than the turn-over mass in Fig. 5, at 
\og{Mc\/ Mq) ~ 5. The arrow in Fig. 5 indicates the mass 
where the onset of signihcant selection effects, and therefore 
of sample incompleteness, is expected to occur, based on 
the selection limit in the (age vs. mass) plane shown in Fig. 
3. The expected effect of the Av ;$ 0.2 mag extinction for 
the clusters with well-determined ages is a marginal shift in 
mass towards higher masses of Alog(Mci/M©) < 0.08. This 
implies that the observed turnover is not a spurious effect 
due to varying extinction. In addition, we do not observe a 
systematic trend between cluster mass and extinction, de¬ 
rived from the SED fitting, which would be expected if the 
higher-mass clusters (and therefore brighter at the same age) 
were located slightly deeper into M82 B, while the fainter 
lower-mass clusters were only observed near the very surface 
of the region. 

Therefore, we conclude that if the initial mass spec¬ 
trum of the M82 B clusters formed in the burst of clus¬ 
ter formation resembled a power-law distribution down to 
the low-mass selection limits, cluster disruption must have 
transformed this distribution on a time-scale of < 10 ® yrs 
into a distribution resembling a log-normal or Gaussian dis¬ 
tribution. 

We now compare the observed mass distribution of the 
clusters formed in the burst with model predictions. These 
predictions are based on the assumption that the clusters 
were formed with a cluster IMF of slope tr = 2, and that 
the distribution has been modihed subsequently by cluster 
disruption. For the cluster disruption we adopt the same 
relation fais = t 4 '®(M/ 10 ^M©)'’' with 7 = 0.62 (Eq. (1)), as 
found in Paper I for four galaxies, and used in the previous 
section. We can then model the mass distribution of the 
observed burst clusters with only one free parameter, i.e. 
the disruption time 14 “ of a cluster with an initial mass of 
10^ M©. 



Figure 5. Mass distribution of the clusters formed in the burst of 
cluster formation, 8.4 < log(t/yr) < 9.4. The shaded histograms 
correspond to the clusters with (log (Age [max]) — log (Age [min]) < 
1.0; the open histograms represent the entire cluster sample in this 
age range. The arrow indicates the mass where the onset of signif¬ 
icant selection effects, and therefore of sample incompleteness, is 
expected to occur, based on the selection limit in Fig. 3. The ex¬ 
pected effect of the Ay < 0.2 mag extinction for the clusters with 
well-determined ages is a shift in mass towards higher masses of 
A log(Mci/M0) < 0.08, which implies that the observed turnover 
is not a spurious effect due to varying extinction. 


The decreasing mass of a cluster with time is given by^ 
dM _ M M 

"dT “ ~ ~4"(M/104Mo)'^ ’ 

so that 

M{Mi,t) = {M^ , (3) 

where M and Mi are, respectively, the present and the initial 
mass of the cluster, both in units of Mq, and B = 7 l 0 '*^/t 4 “. 
Suppose that clusters form in a burst as 

dN{Mi, t) = S{t)M-°‘AMiAt (4) 

over the mass range Mmin < Mi < Mmax, and governed by a 
formation rate S{t) in units of M©““^ yr~^. We assume an 
instantaneous burst of cluster formation, so that S{t) is a 5 
function at age tburst = 10® yr. We adopt Mmax = 3 x 10® M© 
(i.e., logMmax/M© = 6.5, one mass bin beyond the maxi¬ 
mum mass of clusters formed in the burst, thus allowing for 
some disruption effects) and Mmin = 10^ M©, but this lower 
limit is not important since it is below the detection limit of 
^ 10^ Mq for the burst clusters (see Fig. 2). The function 


1 Recent simulations, based on both Monte Carlo realisations 
(e.g., Giersz 2001, his Fig. 2) and A-body modeling (e.g., Porte- 
gies Zwart et al. 2002, their Fig. 2; S.J. Aarseth, priv. comm.), 
show that this is a close approximation to real cluster evolution, 
provided that the initial cluster has reached an equilibrium state. 
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log Mass (Mq) 

Figure 6. Predicted mass distributions in terms of log(d7V/dM) x 
(M/Ntot (equivalent to the logarithm of the number of clusters 
per mass bin) for an instantaneous burst of cluster formation that 
occurred 1 Gyr ago. From top to bottom, the curves represent the 
(normalised) mass distributions for characteristic cluster disrup¬ 
tion times log(^ 4 ^®/yr) = 10.0, 9.5, 9.0, 8.5, 8.0, and 7.5. Notice the 
dependence of both the turnover mass, indicated by the dotted 
line, and the maximum mass on 


S{t) is related to the total number of clusters formed per 
unit time as 

S{t) = ^ (a - 1) j (5) 

where dA 7 ot/dt is the total number of clusters formed per 
year. It is easy to show that the present mass distribution of 
clusters formed in an instantaneous burst of age tburst can 
be written as 

n/T \ 'y \ 

‘ + ''(T03Me) 

for masses above the detection limit and smaller than 
~ Btburst)^^^- We 866 that the initial mass distri¬ 
bution, dN{Mi,t) = has been modified by a fac¬ 

tor that depends on the ratio tburst/tf^, for given values of 
a = 2.0 and 7 = 0.62. 

We have calculated the predicted mass distribution, us¬ 
ing Eq. ( 6 ), for a burst age of 1 Gyr and for several val¬ 
ues of the cluster disruption time-scale tf'®, see Fig. 6 . To 
match the observed peak in the cluster mass distribution, 
a characteristic disruption time-scale log(t 4 '®/yr) between 
7.5 and 8.0 is required, with the most likely value closer to 

log(4“yr) = 8 . 0 . 

We then calculated the expected number of clusters in 
the same age bins as observed (Fig. 5) and normalised the 
distribution to the observed total number of clusters. The 
results are shown in Fig. 7a for the subsample of the 42 clus¬ 
ters with the most accurately determined ages, and in panel 
b for the full sample of 58 clusters. The location of the pre¬ 
dicted peak is very sensitive to the value of tf'^/tburst, which 
follows from Eq. ( 6 ), and is shown in Fig. 6 . We found the 
best fit for a cluster disruption time of 14 '® ~ 3 x 10“^tburst, 
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Figure 7. Comparison between the observed mass distribution of 
the clusters formed in the burst of cluster formation and the pre¬ 
dicted distribution, based on a cluster IMF with slope 0 = 2 and 
cluster disruption with time-scales of Gig = /10‘^Mq)^-^'^ . 

(a) Cluster sample with the most accurate masses and ages, (b) 
Full sample. The predicted number of clusters in each mass bin, 
normalised to the total number of observed clusters, are shown 
by the dotted, full and dashed lines for characteristic disruption 
time-scales log(t 4 '®/yr) = 8.5, 8.0, and 7.5, respectively. 


with an uncertainty of approximately a factor of two. This 
corresponds to fburst = 3 x 10 ^ yr if fburst = 1 Gyr. We con¬ 
clude that the mass distribution of the clusters formed dur¬ 
ing the burst can be explained by an initial cluster IMF with 
a slope of a = 2, an upper limit to their mass of ^ 3 x 10 ® 
Mq and cluster disruption with a time-scale given by Eq. 
( 1 ), ~ 3 x 10 ^ yr if the age of the burst is tburst = 1 x 10 ® 

yr. The uncertainties in both the observed and predicted 
numbers of clusters of a given mass are governed by identi¬ 
cal Poisson-type statistical errors. The uncertainties in the 
mass and age scales, on the other hand, are dominated by 
systematic uncertainties introduced by the adopted stellar 
evolutionary synthesis models, and - to some extent - by the 
specific extinction law adopted, although the relative mass 
and age distributions are robust. 
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6 DISCUSSION 

6.1 The cluster disruption time-scale 

We have derived the characteristic cluster disruption time- 
scale from the mass distribution of the clusters formed dur¬ 
ing the burst. In principle this disruption time-scale could 
also be derived from the the age distributions of the clusters 
formed during the post-burst phase, shown in Fig. 4c. How¬ 
ever the number of observed young clusters is very small, 
only a few per age bin so that the crossing point between 
the (dotted) fading line and the (solid) disruption line is very 
uncertain. However, the crossing point suggested in Fig. 4c 
agrees approximately with that expected for the disruption 
time-scale of ~ 3 x 10^ yr derived from the clusters 
formed during the burst (see the description of the method 
in Paper I). 

In Paper I we have shown that the disruption time-scale 
can also be derived from the mass distribution of clusters in 
a magnitude limited sample if the cluster formation rate 
is constant. This last assumption is certainly not justified 
for M82 B. For a strongly variable cluster formation rate, 
as in M82 B, the mass distribution is heavily affected such 
that the crossing point between the power-law fit to the 
distribution of low mass clusters (the fading line) and the 
power-law fit to the distribution of the massive clusters (the 
disruption line) does no longer represent the value of 
Therefore, we could only use the mass distribution of the 
clusters formed during the burst to derived the disruption 
time-scale. 

An initial comparison among the characteristic cluster 
disruption time-scales found in Paper I reveals large differ¬ 
ences. However, the characteristic cluster disruption time- 
scale for the clusters in region B of its low-mass host galaxy 
M82 is, within the uncertainties, comparable to that in the 
dense centre of the massive grand-design spiral M51. Lamers 
& Portegies Zwart (in prep.) are currently analysing the 
cluster disruption time-scales derived for the galaxies stud¬ 
ied in Paper I and in this paper, using A-body simulations. 
Here, we will therefore simply explore whether this similar¬ 
ity between in the centre of M51 and in M82 B can 
be understood from a comparison of the ambient density in 
both regions. 

In order to estimate the density in M82 B, we deter¬ 
mined its total U-band luminosity, after correcting for a 
nominal extinction throughout the region of Av ~ 0.5 mag, 
Vtot ~ 9.5 from our HST observations (dGOG), so that 
Mv,tot = —18.3, Lv,tot — 1.8 X lO®I/\/, 0 . Spectral synthe¬ 
sis of region B suggests that it has a U-band mass-to-light 
ratio, M/Lv ~ 0.5 — 1.0 (R.W. O’Connell, priv. comm.), 
as expected from the Starburst99 models for a stellar popu¬ 
lation of several 100 Myr, so that the total mass contained 
in the objects providing most of the F-band luminosity is 
M ~ 1.35 X 10® Mq. If we now assume that the volume oc¬ 
cupied by these objects is roughly similar to a sphere with 
a radius determined by the radius of the area in which we 
measured the total flux, we derive an average density for 
M82 B of (p) ~ 2.5M©pc“®, or log(p)(M©pc“®) ~ 0.4. 
Note, however, that we have introduced large uncertainties 
by adopting the above assumptions, in particular because 
we have assumed to have sampled the entire volume of the 
“spherical region” M82 B. This implies that the region has 
been assumed transparent, as opposed to the conclusion in 


dGOG that the M82 B cluster sample contains likely only 
the subset of M82 B clusters on the surface of the “sphere”. 
The implications of our assumptions are therefore that we 
expect to have missed a large number of clusters present in 
the interior of M82 B. For these missed clusters in the inte¬ 
rior we expect the disruption time to be shorter because of 
the higher densities predicted there. 

For the interstellar medium in the centre of M51, 
Lamers & Portegies Zwart (in prep.) derive a mean density 
of (p) ~ O. 6 OM 0 pc“®, based on column density estimates 
by Athanassoula et al. (1987) and similar geometrical argu¬ 
ments as used above for M82 B. Within the large uncertain¬ 
ties involved in such back-of-the-envelope approximations, 
these two estimates of the mean density in the centre of 
M51 and in M82 B are remarkably similar, within an order 
of magnitude, as are their characteristic cluster disruption 
time-scales. Thus, we conclude that, although M82 as such is 
a small, low-mass irregular galaxy, its fossil starburst region 
B has achieved a similarly dense interstellar medium as the 
centre of M51, so that similar cluster disruption time-scales 
are not a priori ruled out. 

6.2 The cluster formation rate 

We have derived the cluster formation history of M82 B from 
the age distribution shown in Fig. 4c. The age distribution 
itself. Fig. 4a already shows a very strong peak around 10® 
years. Part of the steep increase between 8 < log(t/yr) < 9 is 
owing to the use of logarithmic age bins. Therefore we have 
transformed the age distribution into the formation history 
of the observed clusters. Fig 4c. The general decrease with 
age clearly shows the effect of cluster disruption. Disruption 
of clusters formed at a constant formation rate and with a 
mass-dependent disruption time-scale will result in a power- 
law decrease of the observed cluster formation rate (Paper I), 
roughly in agreement with the observations. The deviations 
from the power law, e.g., the peak in the observed formation 
history around 1 Gyr, reflect changes in the real formation 
history, or in the disruption time. 

For magnitude-limited cluster samples with a constant 
cluster formation rate and a constant mass-dependent dis¬ 
ruption time-scale, the disruption time and its dependence 
on mass can be derived from the age distribution. Because 
of the non-constant cluster formation rate and the small 
number of clusters, such an analysis is not possible here. 
Therefore we opted for the reasonable alternative of adopt¬ 
ing the mass dependence of the disruption time-scale, i.e., 
the value of 7 ~ 0.62 in Eq. (1), which seems to be a uni¬ 
versal value found in four galaxies with very different con¬ 
ditions (Paper I). We also assumed that the scaling factor 
t 4 ‘^ of the disruption time-scale is constant. With these two 
assumptions, we derived the ratios between the real clus¬ 
ter formation rates in the pre-burst phase (log(t/yr) > 9.4), 
the burst-phase (8.4 < log(t/yr) < 9.4) and the post-burst 
phase (log(t/yr) < 8.4 of, roughly, 1:2:^. The formation 
rate during the burst may have been higher if the actual 
duration of the burst was shorter than adopted. 

Alternatively, we could have adopted a constant cluster 
formation rate, but wildly variable cluster disruption time- 
scales. In that case the cluster disruption time-scale should 
have been much longer (i.e. much slower disruption) dur¬ 
ing the burst than before the burst, and again much shorter 
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(faster disruption) after the burst. We think that this alter¬ 
native is very unlikely because (a) disruption is a slow pro¬ 
cess that occnrs over an extended period of time, so large 
changes are not to be expected, (b) there is no obvious phys¬ 
ical process that would produce such large changes and (c) 
the alternative of a variable cluster formation rate is much 
more likely, because this effect is observed in several inter¬ 
acting galaxies and in M82 regions A, C and E (dGOG). 
Therefore, we are confident that the changes in the apparent 
formation rates of the observed clusters are due to changes 
in the real cluster formation rate. 


7 SUMMARY AND CONCLUSIONS 

In this paper, we have reanalysed the previously published 
optical and near-infrared HST photometry of the star clus¬ 
ters in M82’s fossil starburst region B (dGOG) to obtain 
improved individual age and mass estimates. We have also 
extended this study to obtain estimates for the cluster for¬ 
mation history and include the importance of cluster disrup¬ 
tion. Our main results and conclusions can be summarised 
as follows: 

(i) Our new age estimates, based on improved fitting 
methods, confirm the peak in the age histogram attributed 
to the last tidal encounter with M81 by dGOG. For the sub¬ 
sample of clusters with well-determined ages we find a peak 
formation epoch at log(tpeak/yr) = 9.04 (tpeak — ITO Gyr), 
with a Gaussian a of Alog(twidth) = 0.27, corresponding to 
a FWHM of A log(t„idth) = 0.64. The corresponding num¬ 
bers for the full sample are log(tpeak/yr) = 8.97 (tpeak — 
0.93 Gyr) for the peak of cluster formation, and for a, 
Alog(twidth) = 0.35 (FWHM, Alog(t„idth) = 0.82). In 
dGOG (see also de Grijs 2001), we concluded that there is 
a strong peak of cluster formation at ~ 650 Myr ago, which 
have formed over a period of ~ 500 Myr, but very few clus¬ 
ters are younger than 300 Myr (log(t/yr) ~ 8.5). Our new 
age estimates date the event triggering the starburst to be 
slightly older. Our estimate of the duration of the peak is an 
upper limit because the observed width of the peak in the 
age histogram may have been broadened by uncertainties in 
the derived cluster ages. 

(ii) The improved mass estimates confirm that the (ini¬ 
tial) masses of the young clusters in M82 B with V < 22.5 
mag are mostly in the range lO'^ — 10 ® M©, with a me¬ 
dian of 10®M© (dGOG). We find that the mean mass of 
our M82 B cluster sample is log(Mci/M©) = 5.03 and 4.88 
for the subsample with well-determined ages and the full 
sample, respectively, corresponding to M© = 1.1 x 10 ® and 
0.8 X 10® M©, respectively. If the initial mass spectrum of the 
clusters formed at the burst epoch was a power law, disrup¬ 
tion effects have transformed it into a broader and flatter 
distribution on time-scales of < 1 Gyr. 

(iii) The (apparent) formation history of the observed 
clusters shows a gradual increase to younger ages, reach¬ 
ing a maximum at the present time. This shows that cluster 
disruption must have removed a large fraction of the older 
clusters. Adopting the expression for the cluster disruption 
time, tdis{M) = t 4 ‘'®(M/ 10 ‘^M©)^ with 7 ~ 0.62, that was 
derived for four galaxies characterized by very different con¬ 
ditions (Paper I), we found that ratios between the real clus¬ 
ter formation rates in the pre-burst phase (log(t/yr) > 9.4), 


the burst-phase (8.4 < log(t/yr) < 9.4) and the post-burst 
phase (log(t/yr) < 8.4) are roughly 1:2:^. The formation 
rate during the burst may have been higher if the actual du¬ 
ration of the burst was shorter than adopted. We see that 
the cluster formation rate in the post-burst phase is much 
smaller than in the pre-burst phase, because the burst has 
consumed a large fraction of the available molecular clouds, 
leaving little material for the cluster formation in the post¬ 
burst phase. 

(iv) The mass distribution of the clusters formed during 
the burst shows a turnover at log(Mci/M©) ~ 5.2. This 
turnover is not due to selection effects, because the mag¬ 
nitude limit of our sample would produce a turnover near 
log(Mci/M©) ~ 4.5. The observed distribution can be ex¬ 
plained by cluster formation with an initial power-law mass 
function of exponent a = 2 up to a maximum cluster mass of 
Afmax = 3 X 10®M©, and cluster disruption given by Eq. (1) 
with the adopted value of 7 = 0.62 with only one free fitting 
parameter: tf’^/tburst — 3 x 10“^. For a burst age of 1 x 10® 
yr, we find that t 4 ‘^ ~ 3 x 10^ years, with an uncertainty of 
a factor of two. 

(v) The time-scale of t 4 ‘^ ~ 30 Myr is much shorter than 
derived in Paper I for any of the SMC (~ 8 Gyr), the so¬ 
lar neighbourhood (~ 1 Gyr), M33 (^ 0.13 Gyr) and even 
shorter than in the inner spiral arms of M51 (~ 40 Myr). The 
characteristic disruption time-scale in M82 B is the shortest 
known in any disc (region of a) galaxy. 
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